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Abstract 

The present study deals with spatial homogeneous and anisotropic lo- 
^ C"| cally rotationally symmetric (LRS) Bianchi-II dark energy model in gen- 

Cl^ eral relativity. The Einstein's field equations have been solved exactly by 

taking into account the proportionality relation between one of the com- 
i^j ponents of shear scalar (<j\) and expansion scalar (#), which, for some 

(3JT^ suitable choices of problem parameters, yields time dependent equation 

of state (EoS) and deceleration parameter (DP), representing a model 
which generates a transition of universe from early decelerating phase to 
present accelerating phase. The physical and geometrical behavior of uni- 
verse have been discussed in detail. 
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1 Introduction 

The discovery of acceleration of the universe stands as a major breakthrough of 
the observational cosmology. The power of observations in cosmology is clear 
from the observations of supernovae of la (SN la) which dramatically changed, 
about a decade ago, the then standard picture of cosmology - of an expanding 
universe evolving under the rules of general relativity such that the expansion 
rate should slow down as cosmic time unfolds. Surveys of cosmologically distant 
SN la (Riess et al. 1998; Permutter et al. 1999) indicated the presence of new, 
unaccounted - for dark energy that opposes the self-attraction of matter and 
causes the expansion of the universe to accelerate. When combined with indirect 
measurements using cosmic microwave background (CMB) anisotropies, cosmic 
shear and studies of galaxy clusters, a cosmological world model has emerged 
that describes the universe at flat, with about 70% of it's energy contained in 
the form of this cosmic dark energy (Seljak et al. 2005). This acceleration 
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is realized with negative pressure and positive energy density that violate the 
strong energy condition. This violation gives a reverse gravitational effect. Due 
to this effect, the universe gets a jerk and the transition from the earlier decel- 
eration phase to the recent acceleration phase take place (Caldwell et al 2002). 
The cause of this sudden transition and the source of accelerated expansion is 
still unknown. In physical cosmology and astronomy, the simplest candidate 
for the DE is the cosmological constant (A), but it needs to be extremely fine- 
tuned to satisfy the current value of the DE density, which is a serious problem. 
Alternatively, to explain the decay of the density, the different forms of dynam- 
ically changing DE with an effective equation of state (EoS), ui = p/p < —1/3, 
were proposed instead of the constant vacuum energy density. Other possible 
forms of DE include quintessence (lo > —1) (Stcinhardt et al. 1999), phantom 
(bj < -1) (Caldwell 2002) etc. While the the possibility lo « -1 is ruled out 
by current cosmological data from SN la (Supernovae Legacy Survey, Gold sam- 
ple of Hubble Space Telescope) (Riess et al. 2006; Astier et al. 2006), CMBR 
(WMAP, BOOMERANG) (Eisentein et al. 2005; MacTavish et al. 2006) and 
large scale structure (Sloan Digital Sky Survey) data (Komatsu et al. 2009), 
the dynamically evolving DE crossing the phantom divide line (PDL) (lo = — 1) 
is mildly favored. Some other limits obtained from observational results com- 
ing from SNe la data (Knop et al 2003) and combination of SNc la data with 
CMBR anisotropy and galaxy clustering statistics (Tegmark et al. 2004) are 
— 1.67 < lo < —0.62 and —1.33 < lo < —0.79, respectively. The latest results 
in 2009, obtained after a combination of cosmological data-sets coming from 
CMB anisotropics, luminosity distances of high red-shift type la supernovae 
and galaxy clustering, constrain the dark energy EoS to —1.44 < lo < —0.92 at 
68% confidence level (Komatsu et al. 2009; Hinshaw et al. 2009) 

Moreover, in recent years Bianchi universes have been gaining an increas- 
ing interest of observational cosmology, since the WMAP data (Hinshaw et al. 
2003, 2007; Jaffe et al. 2005) seem to require an addition to the standard cos- 
mological model with positive cosmological constant that resembles the Bianchi 
morphology (Jaffe et al. 2006a, 2006b; Campanelli et al. 2006, 2007; Hoftuft et 
al. 2009). According to this, the universe should achieve a slightly anisotropic 
special geometry in spite of the inflation, contrary to generic inflationary models 
and that might be indicating a nontrivial isotropization history of universe due 
to the presence of an anisotropic energy source. The Bianchi models isotropize 
at late times even for ordinary matter, and the possible anisotropy of the Bianchi 
metrics necessarily dies away during the inflationary era (Ellis 2006). In fact 
this isotropization of the Bianchi metrics is due to the implicit assumption that 
the DE is isotropic in nature. Therefore, the CMB anisotropy can also be fine 
tuned, since the Bianchi universe anisotropics determine the CMB anisotropics. 
The price of this property of DE is a voilation of null energy condition (NEC) 
since the DE crosses the phantom divide line (PDL), in particular depending 
on the direction. 

The anomalies found in the cosmic microwave background (CMB) and large 
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scale structure observations stimulated a growing interest in anisotropic cosmo- 
logical model of universe. Here we confine ourselves to models of Bianchi-typc 
II. Bianchi type-II space-time has a fundamental role in constructing cosmo- 
logical models suitable for describing the early stages of evolution of universe. 
Asseo and Sol (1987) emphasized the importance of Bianchi type-II Universe. 
Recently Pradhan et al (2011) and Kumar and Akarsu (2011) have dealt with 
Bianchi-II DE models by considering the spatial law of variation of Hubble's pa- 
rameter which yields the constant value of deceleration parameter (DP). Some 
authors (Akarsu and Kilinc 2010a, 2010b; Yadav et al. 2011; Yadav and Yadav 
2011; Kumar and Yadav 2011; Yadav 2011; Adhav et al. 2011 and recently 
Yadav and Saha 2011) have studied DE models with variable EoS parameter. 
In this paper, we presented general relativistic cosmological model with time de- 
pendent DP in LRS Bianchi-II space-time which can be described by isotropic 
and variable EoS parameter. The paper is organized as follows: The metric and 
field equation are presented in section 2. Section 3 deals with the solution of 
field equations and physical behavior of the model. Finally the findings of paper 
are discussed in section 4. 



2 The metric and field equations 

The gravitational field in our case is given by a Bianchi type-II (BII) metric 

ds 2 = —dt 2 + a\(dx\ — x^dx 2 ) 2 + a\dx\ + a 2 dx\, (1) 

with ax, a 2 , &3 being the functions of time only. In what follows, we consider 
the LRS BII model setting a 2 =0,3. 

Given the fact that the dark energy is isotropically distributed, it is enough to 
consider only three Einstein equations (Saha 2011) corresponding to the metric 
(fT|), namely 
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V — a\a\ — \J—g. (3) 
The expressions for expansion and shear for BII metric given by (P) read: 
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and 



! <2l 1 2 3 d 2 1 j rai a 2 zVdil 
cr 1 = a% = al= -1?, ct 2 = z 1 — . (5 

CL\ 3 0-2 o Lai 0*2 0*2 -I 

Let us now define the generalized and directional Hubble parameters. As in 
know, the Hubble parameter was defined by E. Hubble for the FRW model 

ds 2 = -dt 2 + a 2 (dxl + dx\ + ax 2 .), (6) 

as H = Taking into account that y/—g — a 3 it can be defined as H = 



df 

I— L=^EI ; and the directional Hubble parameters as Hi = -^== or Hi = 

1 1 dgu 

2 gu dt ' 

Taking into account that for BII metric (JTJ \f—g = aia 2 a3 = ciia 2 = V and 
gu = af , g 22 = a; 2 a 2 + a 2 , and 333 — a 2 , analogically we define 

rr Ol x§aiai + a 2 a 2 d 2 

1 ~ ^T' 2 ~ ^2 I „2 » ^3 - — U) 
ai T a 2 a 2 

and 

3 K 3 ai a 2 

It should be noted that though for a 2 = ai we have Hi = H2 = H3 as in isotropic 
case, the present definition does not lead to H = (Hi + H2 + i?3)/3. For this 
equality to held, one must set i? 2 = Unfortunately, there is no unique 
definition for directional Hubble parameters. Finally we define the deceleration 
parameter (DP) as 

,= - V -l. (9) 
V 2 

Imposing the proportionality condition, i.e., assuming that the expansion $ 
is proportional to say a\ : 

tfocoj, (10) 
one finds the following relations between the metric functions 

a 2 = a n i, (11) 
with n being some constant. Inserting (lll[) into Q we obtain 

Oi = ^V(an+D j a2 = ^n/(2„+l) ) (12) 

Subtraction of (l2b|) from (P2al) gives 

^_^ + (^) 2 _^i^_4 = 0. (13) 
a 2 ai a 2 a\ a 2 a 2 

Inserting ai and a 2 from (jX3J) into (fl"3"|) we find equation for defining T^: 

y = 2n + 1 y(3-2n)/(l+2n)^ ^ 
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with the solution in quadrature 



, dV = 2n+1 t. (15) 

v / V 4/(2n+l) + C ^2(71-1) 

Eq. ([T5|) imposes some restriction on the choice of n, namely, n > 1. Thus we see 
that the proportionality condition (|10|) in our case does not allow isotropization 
of the initially anisotropic space-time. 

Once V is defined, we can define DP from (J9)) and EoS parameter from 

4(n + l)(2n + 1)1^ - 4(n 2 + 2n)V 2 + (2n + i)2y4/(2n+i) 
4(n 2 + 2n)F 2 - (2n + l)2y4/(2„+i) 
1 4(n+l)(2n+l)W> 

4(n 2 + 2n)V 2 - (2n + l)2y4/(2n+i) 

Thus we see that V plays central role here in defining all physical quantities. In 
what follows we find V from (TTJ} or (fT5"|) for some concrete values of n or C. 



3 Solution of field equations 

One can not solve equation (|15p in general. So, in order to solve the problem 
completely, we have to choose either C or n in such a manner that equation 
(fl5|) be integrable. The easiest way is to set C = in (fl~5|) . In that case one 
dully obtains 



1 / = < 7 ot (2n+l)/(2„-l) 5 C(j 



(2n+ l) 2 



(2n- l)^2(n- 1) 



(2n+l)/(2n-l) 

(17) 



As one sees, in this case V is an increasing function of time, but this solution 
leads to the constant DP. 

Since, we are looking for a model explaining an expanding universe with 
acceleration, we consider the case for nontrivial C, which for a suitable choice 
of n gives the time dependent DP. The motivation for time dependent DP is 
behind the fact that the universe is accelerated expansion at present as observed 
in recent observations of Type la supernova (Riess et al. 1998, 2004; Perlmutter 
et al. 1999; Tonry et al. 2003; Clocchiatti et al. 2006) and CMB anisotropics 
(Bennett et al. 2003; de Bernardis et al. 2000; Hanany et al. 2000) and de- 
celerated expansion in the past. Also, the transition redshift from deceleration 
expansion to accelerated expansion is about 0.5. Now for a Universe which 
was decelerating in past and accelerating at the present time, the DP must 
show signature flipping (see Padmanabhan and Roychowdhury 2003; Amendola 
2003; Riess et al. 2001). So, there is no scope for a constant DP at present 
epoch. So, in general, the DP is not a constant but time variable. 



Thus we consider the Eq. ([15)) with a nontrivial C. For C ^ Eq. (fl5 
allows exact solution only when 4/(2rt + 1) = N, where N is an integer number 
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In this case N can be integer only for n = 1/2 and n = 3/2. Since n > 1 wc 
have only one option left, it is to choose n — 3/2. In this case Eq. (|T5|) reduces 
to 

dV 4* (18) 



which after integration leads 



VvTc 



V = 4t + 20t + 7 

-2--C 



where (3 is the integrating constant and 7 — . 
Inserting equation (|20|) into (fl"2j) . we obtain 



ai = (4* 2 + 2/3* + 7)4 
a 2 = (4* 2 + 2/3* + 7 )i 



(19) 

(20) 
(21) 



The physical parameters such as directional Hubble's parameters (H 1 ,H 2 , H 3 ), 
average Hubble parameter (H), expansion scalar (9) and scale factor (a) are, 
respectively given by 
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(4* 2 + 2/3 + 7) 
a = (4* 2 + 2/3* + 7)3 
The components of shear scalar are given by 

1 6(4* 2 + 2/3* + 7) : 

2 3 4t + /3 
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a;|(4* + /3) £ 3 (4< + /3) 
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The value of DP (g) is found to be 
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Figure 1: Plot of deceleration parameter (q) versus time (t). 
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Figure 2: Plot of EoS parameter (to) versus time (t). 



The sign of q indicates whether the model inflates or not. A positive sign of 
q corresponds to the standard decelerating model whereas the negative sign of 
q indicates indicates inflation. The recent observations of SN la (Riess et al. 
1998, Perlmutter et al. 1999) reveal that the present universe is accelerating 
and the value of DP lies somewhere in the range — 1 < q < 0. Figure 1 depicts 
the variation of DP versus cosmic time as representative case with appropriate 
choice of constants of integration and other physical parameters. 
The energy density of the cosmic fluid (p), EoS parameter (uj) and density 
parameter (fi) are found to be 

a = 21(4t + /3) 2 I (32) 

' 16(4i 2 + 2(3t + 7 ) 2 4(4f 2 + 2/3t + 7) V ; 

= 1 _ 320(4^ + 2ffl + 7 ) 

U 21(8t + 2^) 2 -16(4t 2 + 2^ + 7 ) 1 ' 

63 _ 3(4t 2 + 2^ + 7 ) 

64 16(4i + ,3) 2 1 ' 

From Eq. (j3"3")l follows that at large t when only the quadratic terms stay 
alive, from EoS parameter we find 

320.4< 2 1 320.4 

UJ -> 1 -—75 it = 1 - r = 0, (35) 

21.(8i) 2 - 16.4< 2 (21-1).64 ' V ' 

i.e., under the present assumption the universe is ultimately filled with dust 
only at remote future. 



The initial time of the universe is t = Therefore, at t = 

ff+Vff jh^ foe spatial volume vanishes while all other parameter diverge. 
Thus the derived model starts expanding with big bang singularity at t — 

ff+Vff JI7 w j 1 j c j 1 can sn ift e d to t — by choosing 7 = 0. This singu- 
larity is point type because the directional scale factors ai(t) and 02(f) vanish 
at initial moment. The components of shear scalar vanish at t — > 00. Thus in 
derived model the initial anisotropy dies out at later time. 



Figure 2 depicts the variation EoS parameter (w) versus cosmic time as rep- 
resentative case with appropriate choice of constants of integration and other 
physical parameters. It is shown that the growth of u) takes place with negative 
sign. It should be emphasized that there is a number of models for dark energy 
(quintessence, Chaplygin gas, phantom and many more) and quest for the right 
one is still going on. The main idea for the DE is a negative pressure, so one can 
try with a negative EoS parameter. It should be noted that the quintessence is 
given by a barotropic EoS only with negative parameter. We don't call fluid a 
DE, we just construct DE in analogy with fluid. Figure 3 demonstrates the be- 
havior of density parameter (O) versus cosmic time in the evolution of universe 
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Figure 3: Plot of density parameter (f2) versus time (t). 



as representative case with appropriate choice of constants of integration and 
other physical parameters 



4 Conclusion 

In this paper, we have investigated LRS Bianchi II DE model under the assump- 
tion that & oc <j\ . Under some specific choice of problem parameters the present 
consideration yields the variable DP and EoS parameter. It is to be noted that 
our procedure of solving the field equations are altogether different from what 
Pradhan et al (2011) have adapted. Pradhan et al (2011) have solved the field 
equations by considering the variation law for generalized Hubble's parameter 
which gives the constant value of DP and only the evolution takes place either 
in accelerating or decelerating phase whereas we have considered the propor- 
tionality condition § oc a\ in such a way that gives variable DP which evolves 
from decelerating phase to current accelerating phase (see Fig. 1). Thus the 
present DE model has transition of universe from the early deceleration phase to 
current acceleration phase which is in good agreement with recent observations 
(2006). The model has singular origin and the universe is ultimately filled with 
dust only at remote future. 

The theoretical arguments suggest and observational data show, the universe 
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was anisotropic at the early stage. Here we are dealing not only with the present 
state of the universe, but drawing a picture of the universe from the remote past 
to present day. We use the Bianchi model as one of many models able to describe 
initial anisotropy that dies away as the universe evolves. So though the model 
is anisotropic in the past for small t but it becomes isotropic as t — > oo. In the 
derived model, the EoS parameter (w) is evolving with negative sign which may 
be attributed to the current accelerated expansion of universe. Hence from the 
theoretical perspective, the present model can be a viable model to explain the 
late time acceleration of the universe. In other words, the solution presented 
here can be one of the potential candidates to describe the present universe as 
well as the early universe. 



Acknowledgments 

One of authors (AKY) is thankful to The Institute of Mathematical Science 
(IMSc), Chcnnai, India for providing facility and support where part of this 
work was carried out. Bijan Saha is thankful to joint Romanian-LIT, JINR, 
Dubna Research Project, theme no. 09-6-1060-2005/2013. Finally We would 
like to thank the anonymous referee for his valuable questions which helped us 
to understand the depth of the problem. 

References 

[1] Adhav, K. S., et al: Astrophys. Space Sc. 331, 689 (2011) 

[2] Akarsu, 6., and Kilinc, C. B.: Gen. Relativ. Gravit. 42, 119 (2010a) 

[3] Akarsu, 6., and Kilinc, C. B.: Gen. Relativ. Gravit. 42, 763 (2010b) 

[4] Amcndola, L.: Mon. Not. R. Astron. Soc. 342, 221 (2003) 

[5] Assco, E., and Sol, H: Phys. Rev. D 148, 307 (1987) 

[6] Astier, P., et al: Astron. Astrophys. 447, 31 (2006) 

[7] Bennett, C. L., et al: Astrophys. J. Suppl. Ser. 148, 1 (2003) 

[8] Caldwell, R. R., et al: Phys. Rev. D 73, 023513 (2006) 

[9] Caldwell, R. R.: Phys. Lett. B 545, 23 (2002) 

[10] Campanelli, L., Cea. P. and Tedesco, L.: Phys. Rev. D 97, 131302 (2006) 

[11] Campanelli, L., Cea, P. and Tedesco, L.: Phys. Rev. D 76, 063007 (2007) 

[12] Clocchiatti, A., et al: Astrophys. J. 642, 1 (2006) 

[13] dc-Bcrnadis, P., et al: Nature 404, 955 (2000) 



10 



[14] Eiscntcin, D. J., ct al: Astrophys. J. 633, 560 (2005) 

[15] Ellis, G. F. R.: Gen. Rclat. Gravit. 38, 1003 (2006) 

[16] Hanany, S., et al: Astrophys. J. 545, L5 (2000) 

[17] Hinshaw et al.: Astrophys. J. Suppl. 148, 135 (2003) 

[18] Hinshaw, G., et al: Astrophys. J. Suppl. 180, 225 (2009) 

[19] Hinshaw, G., et al: Astrophys. J. Suppl. 148, 135 (2009) 

[20] Hinshaw, G., et al: Astrophys. J. Suppl. 170, 288 (2007) 

[21] Hoftuft, J., et al: Astrophys. J. 699, 985 (2009) 

[22] Jaffe, J., et al: Astrophys. J. 629, LI (2005) 

[23] Jaffe, J., et al: Astrophys. J. 643, 616 (2006a) 

[24] Jaffe, J., et al: Astron. Astrophys. 460, 393 (2006b) 

[25] Knop, R. K., et al: Astrophys. J. 598, 102 (2003) 

[26] Komatsu, E., et al: Astrophys. J. Suppl. Ser. 180, 330 (2009) 

[27] Kumar, S. and Akarsu, 6.: arXiv: 1110.2408 [gr-qc] (2011) 

[28] Kumar, S. and Yadav, A. K.: Mod. Phys. Lett. A 26, 647 (2011) 

[29] MacTavish, C. J., et al: Astrophys. J. 647, 799 (2006) 

[30] Padmanabhan, T. and Raychowdhury, T.: Mon. Not. R. Astron. Soc. 344, 
823 (2003) 

[31] Permuttcr, S., et al: Astrophys. J. 517, 565 (1999) 

[32] Pradhan, A., Amirhashchi, H. and Jaiswal, R.: Astrophys. Space Sc. 334, 
249 (2011) 

[33] Riess, A. G., et al: Astron. J. 116, 1009 (1998) 
[34] Riess, A. G., et al: Astron. J. 607, 665 (2004) 
[35] Riess, A. G., et al: Astrophys. J. 560, 49 (2001) 
[36] Saha, B.: Cent. European J. Phys. 9, 939 (2011) 
[37] Scljak, et al: Phys. Rev. D 71, 043511 (2005) 

[38] Stcinhardt, P. J., Wang, L. M. and Zlatcv, I.: Phys. Rev. D 59, 123504 ( 
1999) 

[39] Tegmark, M., et al: Phys. Rev. D 69, 103501 (2004) 



11 



[40] Tonry, J. L., et al: Astrophys. J. 594, 1 (2003) 

[41] Yadav, A. K. and Yadav, L.: Int. J. Theor. Phys. 50, 218 (2011) 

[42] Yadav, A. K., Rahaman, F. and Ray, S.: Int. J. Theor. Phys. 50, 871 
(2011) 

[43] Yadav, A. K.: Astrophys. Space Sc. 335, 565 (2011) 

[44] Yadav, A. K. and Saha, B.: Astrophys. Space Sc. 337, 759 (2012) 



12 



